Ab initio structure modeling of complex thin- film oxides: thermodynamical stability 

of TiC/thin-film alumina 



Jochen Rohrer, 1 '^ Carlo Ruberto, 1,2 and Per Hyldgaard 1 

1 BioNano Systems Laboratory, Department of Microtechnology, 
MC2, Chalmers University of Technology, SE-412 96 Gothenburg 
e Materials and Surface Theory Group, Department of Applied Physics, 
Chalmers University of Technology, SE-412 96 Gothenburg 
(Dated: March 6, 2009) 

We present an efficient and general method to identify promising candidate configurations for 
thin-film oxides and to determine structural characteristics of (metastable) thin-film structures us- 
ing ab initio calculations. At the heart of this method is the complexity of the oxide bulk struc- 
ture, from which a large number of thin films with structural building blocks, that is motifs, from 
metastable bulk oxide systems can be extracted. These span a coarse but well-defined network of 
initial configurations for which density functional theory (DFT) calculations predict and implement 
dramatic atomic relaxations in the corresponding, resulting thin-film candidates. The network of 
thin-film candidates (for various film thicknesses and stoichiometries) can be ordered according to 
their variation in ab initio total energy or in ab initio equilibrium Gibbs free energy. Analysis of 
the relaxed atomic structures for the most favored structures gives insight into the nature of stable 
and metastable thin-film oxides. We investigate ultrathin alumina nucleated on TiC as a model 
system to illustrate this method. The stable a- and metastable K-AI2O3 bulk structures lead to an 
alumina-film candidate-space that consists of 38 configurations for a given film thickness, includ- 
ing three different stoichiometries. We identify the stoichiometries that are relevant in equilibrium 
with an O environment from ab initio thermodynamics calculations of the relaxed configurations. 
These relevant stoichiometries are ALi n _406n and ALt n _206n (only in equilibrium at extremely low 
O chemical potentials), with n = 2, 3, 4 identifying the number of oxygen layers. The films with 
AUnOen stoichiometry are not stable for any allowed value of the O chemical potential. Our analysis 
of the atomic structure shows that the favorable structural motifs of the relaxed films heavily differ 
from those in the bulk. In particular the number of tetrahedrally coordinated Al ions is much higher 
in the films and the corresponding tetrahedra are oriented differently than in the bulk. This finding 
of additional or novel favorable motifs documents that the method is capable of catching thin-film 
candidates with a structural nature that is not explicitly included in the network of initial thin-film 
configurations. Our analysis also shows that the thermodynamically most stable TiC/Al4„_406n 
systems decay into a partly decoupled TiC/0/Al4„_4 06n-6 system, with only a weak binding of 
the Al4 n _406n-6 film on the TiC/O substrate. 



PACS numbers: 68.55.-a, 68.47.Gh, 68.35.-p, 64.75.St 

I. INTRODUCTION 

Understanding the atomic and electronic structure of 
thin-film oxides is of significant industrial and fundamen- 
tal importance and a huge challenge at the same time. 
Bulk oxides are characterized by a strong ionicity, which 
often results into a tendency for a high structural flexibil- 
ity and an organization in a large number of different sta- 
ble and metastable phases. Prominent examples can be 
found among aluminum oxidesji titanium oxides^ vana- 
dium oxides^ or hafnium oxides^ For an ultra-thin film, 
the structural variety of the oxide can be even larger 
The mainly insulating character of oxides makes accu- 
rate experimental atomic and electronic structure de- 
terminations difficult, since high-resolution techniques 
(low energy electron diffraction (LEED), 5 scanning tun- 
neling microscopy (STM), 6 transmission electron mi- 
croscopy (TEM)j£ scanning electron microscopy (SEM), 8 
. . . ) mainly use charged particles. Theory assisted meth- 
ods, such as density functional theory (DFT) calcula- 
tions, are of high complementary value. However, when 



modeling thin films that are adsorbed on a substrate, 
relatively large surface unit cells are often needed. As 
a consequence, an enormous number of possible atomic 
configurations for the film arises, and a structure determi- 
nation by straightforward energy calculations of all pos- 
sible candidates becomes computationally intractable. 

The nucleation of alumina on TiC provides an illustra- 
tion of the complexity and importance of predicting and 
understanding atomic structure in oxides, ultra-thin ox- 
ide films and their interfaces. Multilayers of TiC/alumina 
are highly relevant for industrial application as wear- 
resistant coatings on cemented-carbide cutting tools^ 
They are commonly fabricated by chemical vapor de- 
position (CVD). Typically, the a-A^Oa (stable in the 
bulk) and K-AI2O3 (metastable in the bulk) phases are 
obtained with relative orientations are cv(0001)| TiC(lll) 
and /«{001}||TiC(lll) J£ However, these ordered struc- 
tures only arise when the alumina possesses a consider- 
able thickness. The nucleation of alumina on TiC in- 
volves the formation of ultra-thin alumina films. In- 
sight into the detailed, atomic configuration in the 



ultra-thin films is essential because their structure may 
strongly influence the subsequent growth^ A complete 
search through all possible atomic thin-film configura- 
tions by total energy calculations is, however, extremely 
difficult^ 

Of course, ab initio molecular dynamics (MD)i^ is 
a powerful tool that can generally sample typical, and 
hence relevant, structural configurations (for a given film 
thickness and stoichiometry) by simply following the 
atomic dynamics at some elevated temperature. The use 
of density functional theory (DFT) to evaluate the forces 
on the constituent atoms makes this a highly accurate, 
but also costly, method. For strongly ionic materials like 
(thin-film) alumina drastic charge rearrangements occur 
with the motion of every single atom. As this calls for 
very small time steps, formidable ab initio MD simula- 
tion times can be expected. Moreover, for a study of 
thin-film nucleation it is imperative to explore possible 
candidate structures for a range of different thicknesses 
and stoichiometric compositions. It is clear that some 
alternative, accelerated search method is desirable. 

In this paper, we suggest an ab initio method to search 
for the structural elements of thermodynamically sta- 
ble and metastable thin-film oxide configurations nu- 
cleated on metallic substrates. The ab initio method 
avoids costly molecular-dynamics simulations but invokes 
atomic relaxations specified by underlying DFT calcula- 
tions of the atomic forces. The approach makes use of 
the complexity of oxide bulk structures and their ionicity. 
The method is illustrated for the TiC/thin-film alumina 
system that we investigate as an example. We start with 
a coarse sampling of the thin-film configuration space by 
relaxing all possible films that consist of partial alumina 
bulk structures. Due to their stability in the bulk they 
can be assumed as promising initial configurations that 
contain structural elements, that is motifs, relevant for 
the actual stable and metastable thin-film structures. To 
further trigger relaxations towards other configurations 
that do not possess a partial bulk structure, we slightly 
distort the unrelaxed candidates. The relaxed atomic 
structures that are energetically favorable and that differ 
essentially from the partial bulk structures can be used 
to design new candidates with a similar structure. The 
procedure could then be iterated until self-consistency is 
reached. For thin-film configurations that differ in their 
stoichiometry, the energy criterion is replaced by the cri- 
terion of lowest Gibbs free energy. The latter we calculate 
from ab initio thermodynamics! 14 ' 15 i 16 

The paper is organized as follows: Section [II] summa- 
rizes the properties of alumina and TiC that are rele- 
vant for TiC/thin-film alumina. In Sec. IIIII we derive 
all TiC/thin-film alumina initial configurations that are 
consistent with the bulk structure of the respective ma- 
terials. The details concerned with the computation of 
total and Gibbs free energies are discussed in Sec. [TV] In 
Sec. El we present our results on the energetics and ther- 
modynamical stability of thin-film alumina. An analysis 
of the atomic structure of relaxed films is given in Sec. I VII 
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Color coding for stacking sites 

FIG. 1: Bulk structures of a- (left) and K-AI2O3 (right) within 
orthorhombic unit cells. The top panels show side views 
along [100] . The bottom panels define the atomic site labeling 
within each (001) atomic layer. 
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TABLE I: Mapping of stacking and site labels (as defined in 
Fig. [TJ under mirror transformations. A reflection about the 
x z- plane leads for example to a relabeling of cp — > c 7 . 



In Sec. lVII| we discuss our results and Sec. lVIIll contains 
our conclusions. 



II. MATERIALS BACKGROUND 

We first summarize the bulk and surface properties 
of alumina and TiC that are relevant for our thin-film 
modeling method. 



A. Stable and metastable AI2O3 bulk structures 

Figure Q] shows a schematics of the bulk structures of 
a (trigonal unit cell, space group i?3c) and k alumina 
(orthorhombic unit cell, space group Pna2i) i 17 i 18 Along 
the a[0001] and k[001] directions, both alumina phases 
are composed of alternating O and Al layers, the latter 
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splitting up into two sublayers. In a-A^Os, all Al ions 
are octahedrally (O) coordinated. In k-A1 2 3 , the co- 
ordination alternates. In every second layer all Al ions 
have octahedral coordination. In the other layers 50 % 
of the Al ions is octahedrally and 50 % tetrahedrally (T) 
coordinated. All tetrahedra point in the [001] direction. 

To facilitate a parallel treatment of a- and K-AI2O3, we 
choose an orthorhombic unit cell for the representation 
of both alumina phases, so that a[0001] hcx 4=> a[001] ortho . 
The associated calculated lattice parameters are a = 
4.798 (4.875) A, b = 8.311 (8.378) A, and c = 13.149 
(9.018) A for a (k)^S. which is in good agreement with 
experimental data. 20 ! 21 

The atomic structures of a- and K-AI2O3 can be de- 
scribed as follows. We denote the stacking sites of full O 
layers by capital letters. For the Al layers small letters 
with subscript (Arabic numerals for a, Greek letters for 
k) are used. The subscript relates each two Al sites per 
unit cell, see Fig. [1] for a detailed definition of each label. 
For k this notation is identical to the one introduced in 
Ref. [HI, whereas for a it is a slightly modified version of 
the one of Ref. where Greek superscripts are used for 
the labeling of Al vacancies. With this notation, the bulk 
stackings ar o 18 i 22 

a[0001] : Ac3C2BciC3Ac2CiBc3C2Acic 3 Bc 2 ci 
k[001] : Ab 1 cpBc a c 1 Acf3b 1 Cb a bi3 . (1) 

Reflections about the zz-plane [<=> (010)] or yz-plane 
[<;=> (100)] are symmetries of the bulk. Although the 
structure is not invariant under these transformations, 
the transformed structures are equivalent to the non- 
transformed one. The effects of the mirror transforma- 
tions on individual stacking sites are listed in Table [U 
Reflection about the xz-plane corresponds to interchang- 
ing <-y 7 and 2^3, reflection about the j/z-plane to 
C <-> B and c «-> b. 



B. TiC(lll) surface and reactivity 

Bulk TiC possesses NaCl structure with a theoretical 
lattice parameter— a = 4.332 A (in good agreement with 
the experimental valued a exp = 4.33 A). Along the [111] 
direction, it is thus composed of close-packed alternating 
Ti and C layers. The stacking sequence of one repeat 
unit is ABC ABC. 

We only consider Ti-terminated TiC(lll) surfaces. 
This choice is motivated by the stronger binding of Ti to 
the C-terminated surface compared to the binding of C 
to the Ti-terminated surfaced Furthermore, there is ex- 
perimental evidence for a preferred Ti termination upon 
annealing^ 

On Ti-terminated TiC(lll), atomic O adsorbs much 
more strongly than atomic Al (about three times as 
strong)^ We therefore identify the first alumina layer 
above the TiC/Al203 interface plane as an O layer. Ac- 
cording to Refs.[23landl26l. both single O atoms and a full 



O monolayer prefer adsorption in the fee site. By defin- 
ing the TiC stacking such that the fee site on its (111) 
surface is labeled by an A stacking letter, the hep site 
by B, and the top site by C, the position of the first O 
layer is therefore fixed to A stacking. For the monolayer, 
our calculated Ti-0 layer separation along TiC [111] is 
rf Ti _o = 0.89 A. 



III. THIN-FILM IDENTIFICATION METHOD 

We seek a characterization of alumina nucleation on 
TiC and identify promising thin-film alumina structure- 
candidates for different oxide layer thicknesses and sto- 
ichiometrics. This section describes in detail the pro- 
posed method for finding promising initial configurations 
for TiC/thin-film alumina candidates. First, we derive 
all the TiC/alumina interface configurations that are 
consistent with the respective bulk structures and that 
take into account the adsorption properties of TiC(lll). 
Then, we obtain all the initial thin-film configurations 
that consist of partial bulk alumina by truncating these 
interface sequences. 

A. TiC/alumina interface structures 

Table H] lists all the TiC/a-Al 2 3 and TiC/K-Al 2 3 
interfaces that are conform with the respective bulk 
structures and that start with an O layer in fee (A) site 
on the Ti-terminated TiC (111) surface. These stacking 
sequences are found as follows: 

We observe that any of the O layers in the listing (fTJ) 
can be chosen as the initial alumina layer. This layer 
must be translated to an A site, which can be achieved 
by cyclic permutations. For the C sites to be translated 
to A sites we need one cyclic permutation, for B sites 
to be translated to A sites we need two. All other sites 
are relabeled accordingly. For example, for C — > A, we 
have A — > B and B — > C. For the Al positions the 
corresponding relabeling has to be performed, keeping 
the subscripts [a (1), . . .] fixed. 

Next, we note that the fixed stacking sequence of the 
TiC substrate breaks the symmetry associated with a 
reflection about the yz-plane. Hence, for each alumina 
sequence, we need to consider an additional one, which 
is obtained by interchanging B{b) «-> C(c). 

Finally, we exploit that reflection about the rrz-plane 
is still a symmetry of TiC / alumina since the TiC is com- 
posed of fully occupied layers. Hence, alumina sequences 
that are related by <-> 7 (2 «-> 3) are equivalent. 

For a-Al 2 03, all O layers are equivalent. Therefore it 
is sufficient to focus on the first O layer, which is already 
in A stacking. Also, a[001] a[001] and thus only the 
symmetry breaking associated with the reflection about 
the yz-plane needs to be considered. As a result, only 
two possible interfacial configurations have to be taken 
into account (see Table HT1 left column). 
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TiC/g[0001] TiC/ k [001] TiC//t[OQl] 

TiC/ AC3C2BC1C3AC2C1BC3C2AC1C3BC2C1 TiC/ Ab^cpBcaC-fAcpb^Cbabp TiC/ Ab^baCb^cgAc^CaBcgb^ 

TiC/Ab2b 3 Cbib2Ab 3 biCb2b 3 Abib2Cb s b 1 TiC/ Ab a b 7 Cbga^Ba a agCa~fbi3 TiC/ Abga^CagaaBa^bgCb^bc 

TiC/ Acgb 1 Cb a bgAb 1 cgBc a c 1 TiC / Ac^CcBcgb-yAbgbaCb^cg 

TiC/ Ac a cgBc 1 agCa a a, 1 Bagc 1 TiC/ ' Ac 1 agBa 1 a a Cagc 1 Bcgc a 



TABLE II: TiC( 111) /alumina interface configurations that respect the bulk structure of a- and K-AI2O3 and start with an O 
layer in fee site on the Ti-terminated TiC(lll). The TiC stacking of the surface region is defined as . . . ABC ABC . 



For K-AI2O3, only every second O layer is equivalent 
and k[001] is not equivalent to «[001]. Therefore, we need 
to consider both directions, any two consecutive bulk O 
layers, and the effect of the symmetry breaking. This 
results in four different configurations for each direction 
(see Table HH middle and right columns). 



B. TiC/thin-film alumina candidate structures 

We obtain the network of initial thin-film alumina con- 
figurations in three steps. In the first step, we truncate 
the TiC/alumina interface sequences in Table HT1 after a 
full Al layer. The number of O layers n defines the thick- 
ness of the film. In a second step, the resulting configu- 
rations are distorted by placing the Al sublayers into one 
and the same plane, exactly in between the two neighbor- 
ing O layers. In the third step, we vary the stoichiometry 
by removing Al ions from the surface in accordance with 
the bulk space group, i.e., only Al pairs that belong to 
the same stacking label are removed. 

In this way, for each thickness we generate three stoi- 
chiometry classes: Ai4„06n, Ai4„_206n, and Ai4 n _406n, 
corresponding to the removal of zero, one, and two Al 
pairs, respectively. 

For the a-Al 4 „_2 06n films, we only consider the sur- 
face Al pair that has no direct neighbor in the layer be- 
low. With this choice, we minimize the electrostatic en- 
ergy. We have confirmed the quality of this choice by 
comparing the total energy for some configurations with 
different choices of the surface Al pair. 

For the K-Al 4 „_2 06n films, there is no such simple ar- 
gument and we choose to allow for both possible Al pairs. 
As a result, the number of K-Ai4„_206n configurations 
is twice as large as the number of K-Al 4 „_ 4 06n or k- 
Ai4„06n configurations. 

In Tables IIIIfFVl all candidates found by the described 
procedure are listed for films with thicknesses n = 2, 3, 
and 4 (for thermodynamical reasons, see Sec. IV Bl no 
AI16O24 configurations are considered). 



IV. AB INITIO METHOD 
A. Total energies and atomic relaxations 

All calculations are performed with the DFT plane- 
wave code dacapo2£ using ultra-soft pseudopotentials 28 
and the PW91 exchange-correlation 29 functional. 

We use a supercell approach and model the TiC /thin- 
film alumina by slab geometry. The basal plane dimen- 
sions of the supercell are chosen to fit the 3x2 TiC(lll) 
surface (5.306 x 9.190 A 2 ) and the height is fixed to 30 A, 
ensuring a vacuum thickness of at least 13 A. 

The TiC is modeled by four atomic bilayers (with six Ti 
and six C atoms per bilayer). The alumina films contain 
six O atoms per O layer and a varying number of Al 
atoms, depending on the film stoichiometry. In total, the 
slabs contain between 64 (AI4O12 films) and 86 atoms 
(Ali 4 24 films). 

We use a 400 eV plane- wave cutoff and a 4x2x1 
Monkhorst-Pack 3 -^ fc-point sampling. Electrostatic ef- 
fects arising from the charge asymmetry in the slab are 
corrected for by a dipole correction. The atomic re- 
laxations are performed until all interatomic forces are 
smaller than 0.05 eV/A. This choice has proven a good 
accuracy at acceptable CPU times for a- and K-AI2O3 
surfaces^ and for TiC/alumina interface calculations. 31 
The presented DFT calculations amount to a total of one 
million CPU hours on modern supercomputing facilities. 

B. Equilibrium Thermodynamics 

At non-zero temperature T and pressure p the stability 
of any system is governed by the Gibbs free energy G. 

Surfaces. Figure HJa) shows a typical atomic 
setup used for calculating surface energies using slab 
geometry* 3 ^ For alumina, the average surface Gibbs free 
energy cr av of the pair of alumina surfaces represented by 
the slab is defined by 

o- av = i (er + + er~) = i (G 3lab - n A1 ^ A1 - UoHo) . (2) 

Here a + and a~ correspond to the two surface energies 
associated with each side of the slab (cr av = a + = a~ = a 
for a symmetric slab) , G slab is the Gibbs free energy of the 
slab that contains n Al Al and n Q O atoms, and fi AJ and 
jj, are the chemical potentials of Al and O respectively. 
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FIG. 2: Schematics of atomic setup for calculations of (a) (ox- 
ide) surfaces (b) (metal/oxide) interfaces and (c) thin-film 
oxide on a metal substrate. The arrows point to the regions 
in the slabs where the bulk (g bulk ), surface (a), interface (7) 
or thin-film (g fllm ) contributions to the Gibbs free energy are 
located. We emphasize that, in general g^lda / SoxkL- 



The first two contributions can be expressed by 
^(7 + a) 



TiC 



(5) 



Because we always consider the same TiC slab (always 
having identical surface energies) it is convenient to ex- 
press n TiC g TiC through the calculated Gibbs free energy 
of a TiC slab 



G-riC — ^TiC.f?TiC "T" A(7 1 



TiC 



Aa 



TiC 



(6) 



The stability of the thin-film system is therefore conve- 
niently described by the Gibbs free energy difference 



r 



A( 7 + cr 

G T iC/alum: 



' TiC 



2a+ c ) 



TiC/ 

^aiMai 



(7) 



In equilibrium with an O environment, the stoichio- 
metrically weighted sum of the Al and O chemical poten- 
tials must equal the Gibbs free energy per stoichiometric 
unit of alumina 5ai 2 o 3 , 



2/^Al + 3^t — 3ai 2 o 3 



so that CT av can be rewritten as a function of the chem- 
ical potential only. 

Interfaces. Figure [2Kb) shows a schematics of 
of a typical slab geometry used to calculate interface 
energies^ For metal/alumina interfaces, the stability is 
determined by the interface Gibbs free energy 7. Us- 
ing the slab geometry of Fig. H^b), the average interface 
energy is calculated as 




riifii 



Aa+ 

^ u metal 



c t a 1 



„ ffAl 2 3 



Mo ^mctal ~ ^°metal 



(4) 



Here 7 + and 7~ are the interface energies corresponding 
to the two interfaces in Fig. Ufa), G slab is the Gibbs free 
energy of the total slab, and rii and [ii are the number and 
chemical potentials of the different atomic species (i £ 
{Al, O, or metal}). In the last equality, we have assumed 
equilibrium with an O environment and used Eq. [3] to 
express 7 av as a function of the O chemical potential only. 
For a symmetric slab, we have 7 av = j + = 7~ =7. 

TiC '/thin-film alumina. Figure [He) sketches the 
atomic setup used for the present thin-film calculations. 
The stability of TiC/thin-film alumina is governed by 
three contributions: (i) the surface stability of the thin 
film (a) , (ii) the stability of the interfacial configuration 
(7), and (iii) the internal stability of the thin film itself 
(9 ! 1 



film 

'AloOa 



Note that for practical reasons we do not normalize T by 
the area of the basal plane of the unit cell, since our unit 
cells all have the same basal plane. 

For the third contribution that governs the stability of 
thin films, g A1 m , we note that in general the Gibbs free 



(3) energy per stoichiometric unit may and is expected to 



differ from that in the bulk, that is, ^ 9a\ 2 o 3 - We 



therefore introduce a parameter 6, which measures the 
difference between the Gibbs free energy of one stoichio- 
metric unit of alumina in a bulk environment and in the 
film and rewrite the equilibrium condition as 



2fi Al + 3/Lt c 



+ 5 . 



(8) 



The stability-determining quantity for TiC/thin-film alu- 
mina configurations can be reformulated from Eqs.[7]and 

m 



r = g 



TiC/Al„0„ 



G T iC (ffAl 2 3 



8) 



-(n c 



2 "Al) MO 



(9) 



The limits of the physically allowed range of the chemical 
potentials are defined by Al condensation into fee Al and 
O condensation into O2, i.e. /i A1 < g fcc _ A1 and < 
h/J,o 2 , where g fcc _ A i and ^ D2 are the Gibbs free energy per 
stoichiometric unit of fee Al and the chemical potential of 
O2 , respectively. Combining both inequalities and Eq. [5] 
yields 



2 (<?ai 2 o 3 



S - 2g foc . A 



1 

< Mo < 2^0. 



(10) 



Although an exact value of S cannot be calculated, we 
can estimate S by calculating energy differences between 
films that differ by an integer number of stoichiometric 
units, 



E 



Al Tl _ 4 0„ 



2c 



A1 2 3 



)/2 , (11) 



We find that S nm is 0.4 eV and 0.7 eV when compar- 
ing three and two and four and three layer thick films 
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for Al4„_406n stoichiometry. For Ai4„_206n stoichiom- 
etry the corresponding values are 0.3 eV and 1.2 eV. For 
Al4„06n stoichiometry we have only considered three and 
two layer thick films for which we find S nm = 1.1 eV. 

The fact that the largest values of S nm are found when 
calculating the energy differences for the thickest consid- 
ered films is counterintuitive. We would expect that the 
difference in Gibbs free energy per stoichiometric units 
converges towards that of the bulk once the film is thick 
enough. This shows the difficulties in determining the 
Gibbs free energy of a thin film properly. The higher val- 
ues for thicker films may be due to completely different 
surfaces of the respective films and thus due to surface 
energies. 

In the following, we disregard the fact of a non-zero 
value and the stoichiometry and thickness dependence 
of 5, that is, we put 8 = 0. We have checked that the 
resulting uncertainty in T, although certainly not negli- 
gible, does not change our qualitative results as long as 
the temperatures are not too high (below 1300 K). 



V. RESULTS I: ENERGETICAL AND 
THERMODYNAMIC AL STABILITY 

In Tables ITOHVl we list all two, three, and four O layer- 
sthick (initial) alumina-film configurations that are con- 
sistent with the bulk a- or K-AI2O3 structure. The con- 
figurations are grouped according to their stoichiometry, 
film thickness, and the phase and orientation of the alu- 
mina bulk structures from which they are derived. 

For each configuration, we list the stacking sequence 
and the coordination of the Al ions of the alumina film 
before relaxation, together with the calculated total en- 
ergies after relaxation. These energies (E IBl ) are given 
relative to the energy Eq of the energetically lowest ly- 
ing configuration of the same thickness and stoichiometry 
class, E tei = E - E . 



C. Ab initio Equilibrium Thermodynamics 



A. Energetics and metastability 



For the calculation of the Gibbs free energies and chem- 
ical potentials involved in Eq. [9] we use the method de- 
scribed in Ref. 1161 . Here, we only summarize their argu- 
ments and give the computational prescription. 

The Gibbs free energy of solid material, that is of the 
bulk or a slab, is essentially independent of the pressure. 
The temperature dependence is considerably stronger 
and larger in absolute value. Based on the calculated 
vibrational surface Gibbs free energy for Ru02r^ we es- 
timate the vibrational Gibbs free energy per cell for alu- 
mina at T = 1000 K to T vib ~ 1 - 2 eV. Compared to the 
uncertainty in T due the uncertainty in 5, the vibrational 
contributions can therefore savely be neglected.— 

We consequently choose to approximate T by 



r 



Et[ 



C/A1„0„ 



(n - 2 n Ai) Mo , 



(12) 



where E TiC/Aln0m > E TiC , and e A i 2 o 3 are DFT total energies 
of a TiC/Al„O m slab, an isolated (clean) TiC slab and 
one stoichiometric unit of bulk alumina respectively. 

To estimate the partial O2 pressure and the tempera- 
ture that correspond to different values of fJL , we adopt 
the ideal gas approximation. This approximation allows 
us to rewrite fj, as 



Mo(T,p) 



1 



e 02 T + < 5 Mo 2 C r ,p ) 



k B T In 



P_ 

Pa 



.(13) 



Here e^ T is the DFT total energy of the O2 molecule and 
5/j,(T,po) is related to the entropy S and enthalpy H at 
a fixed pressure po, see Ref. [l6| for details. To calculate 
5fj,(T,po), we use the values of S and H for different 
temperatures at standard pressure po = 1 atm that are 
tabulated in Ref. 13a 



For each thickness and stoichiometry class we use the 
total energy of the energetically most-favorable configu- 
ration (the candidate for the stable thin film structure) to 
define a zero-point of relative energy differences, E rol = 0. 
We stress, however, that bulk alumina also exists in a 
number of metastable phases. We may therefore also ex- 
pect metastable configurations among the thin films. 

In our calculation for bulk alumina, we find that the 
DFT energy difference between the k and a phases is 
A QK ~ 0.7 eV/Al 2 03. We use this quantity as an indica- 
tive measure of the metastability of the alumina films 
and define -E mota = 2L Q A aK , where L Q is the number of 
O layers in the alumina film. For stoichiometric films, 
2L Q is equal to the number of stoichiometric AI2O3 units 
in the film. For non-stoichiometric films, it will give 
an approximate measure of the number of stoichiometric 
units. We then consider configurations with E tel > E meta 
as unstable and configurations with < E rcl < E meta as 
potentially metastable. Whether or not they are truly 
metastable cannot, however, be inferred from our calcu- 
lations. Such an analysis would require a calculation of 
the nature of vibrational excitations and is much beyond 
the present investigation. 

For the two, three, and four O layers thick films the 
criterion is E meta = 0.28 eV, 0.42 eV, and 0.56 eV, re- 
spectively. In Tables IIIIfTVl the configurations that are 
stable or potentially metastable are underlined. For the 
Al4n-406n stoichiometry class, there are at most two po- 
tentially metastable configurations for each alumina film 
thickness. For Ai4„06ri, no potentially metastable con- 
figurations are found. Among the Al4„_206n configura- 
tions, metastable films are more common. However, their 
number decreases as the film thickness increases. 
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TABLE III: Stacking sequence and Al coordination [O for octahedral, T for tetrahedral, with the arrow indicating the direction 
in which each tetrahedron vertex is pointing: towards the film surface (f) or towards the TiC/film interface (J.)] of unrelaxed 
alumina films with Al4 n _406n stoichiometry and their total energies E Ici after relaxation (given relative to the structure with 
lowest total energy for each film thickness). The configurations are grouped together according to the phase and orientation of 
the alumina bulk structures from which they are derived (left column). Configurations that differ only by a rotation of 180° 
around TiC[lll] are organized into subgroups separated by larger whitespace. In general, the unrelaxed and relaxed atomic 
structures differ considerably. The stable and potentially metastable (see text for details) configurations are underlined. The 
ab initio study and comparison permit us to make the following set of observations: (i) The unrelaxed configurations with 
an AC stacking in the first two O layers yield relaxed structures that are in general more favorable than those obtained from 
configurations in which the stacking sequence has been rotated by 180° around TiC[flf] {AB O stacking); (ii) While for the 
AI4O12 films two different unrelaxed structures lead to the stable configuration, for the other two film thicknesses only one 
structure leads to the stable configuration; (iii) In general, the stable configurations are obtained from TiC[lll]/ft[001] initial 
structures; (iv) The Q-type films lead to neither stable nor metastable configurations; and (v) While the stable AI4O12 and 
AlgOi8 films are both obtained from the same unrelaxed interface sequence (same line), the stable AI12O24 film derives from 
another interface sequence. 




FIG. 3: Thermodynamic stability of thin-film alumina with stoichiometrically different compositions and different thicknesses 
on the TiC(lll) substrate in equilibrium with an O2 environment. The left panel shows the Gibbs free energy differences V 
fEq. II2[) per unit cell of TiC/thin-film alumina as a function of the O chemical potential A^io = fJ-o ~ | e o 2 (Eq.ll3p for all three 
considered thicknesses and stoichiometric compositions. The left end of each line is defined by the physically allowed range 
(fcc-Al condensation) of the O chemical potential (Eq. I10p . For all thicknesses the alumina films with Al4„_406„ stoichiometry 
(solid lines) are stable at medium to high O chemical potential, whereas films with Al4„_206n stoichiometry (dashed lines) are 
stable at low O chemical potential (A/Xo < —3.5 eV). The alumina films with Al4„06n stoichiometry (dashed-dotted lines) are 
not stable at any allowed value of the O chemical potential. The right panel shows the O chemical potential A/j, as a function 
of partial O2 pressure for three different temperatures. We find that a value of A/i < —3.5 eV at a temperature T = 1300 K 
corresponds to an O2 pressure of Po 2 ~ 10 -15 bar. 
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"The stacking direction ([001] or [001]) cannot be inferred at the 
considered thickness. 



TABLE IV: Stacking sequence and Al coordination of unrelaxed alumina films with Al4„_20Bn stoichiometry and their relative 
total-energy differences J5 ro i after relaxation. Notation and grouping are the same as in Tab. Mil Configurations that differ 
only in their surface Al ion are grouped together and separated by larger whitespace. The coordination given for the surface 
Al ion is the one that it would have in the bulk. The ab initio study and comparison permit us to make the following set of 
observations: (i) Although the stable films are generally of k type, a-type films are competitive, at least for the thinnest films; 
(ii) For the thinner films, both «[001] and «:[00l] orientations yield stable and metastable configurations, while for the thicker 
films, only k[001] leads to (meta-)stable configurations; (iii) The general trend in stability with respect to the O stacking is 
the same as for the Al4„_406n films (AC more favorable than AB) but with exceptions, in particular, the stable AlioOis and 
AI14O24 configurations originate from structures with AB stacking in the first two O layers. 



B. Thermodynamical stability 

Figure [3] shows our calculated values of T for the ener- 
getically most favorable configurations of each considered 
thickness and stoichiometry class. Corresponding values 
for the potentially metastable configurations can be ob- 
tained by adding the relative energies of Tables IlllrfVl 

In the physically allowed range of the O chemical po- 
tential /jL (Eq. ITUj) . we find that the stable film belongs 
to either the Al4 n _206n or the Ai4 n _406n stoichiome- 
try class, independent of the thickness. The Ai4„_406n 
films are stabilized for A/i > —3.5 eV, whereas for 
A/x < —3.5 eV, the Ai4„_206n films are stable. These 
findings apply when the films are in equilibrium with an 
O2 atmosphere. 

In the right panel of Fig. [3] we show the relation be- 
tween the O chemical potential and the O2 pressure at 
several temperatures. We find that an O chemical po- 
tential of A/i Q > —3.5 eV can only be reached for rela- 
tively high temperatures and extremely low O2 pressures 
(T - 1300 K, po 2 ~ 10~ 15 bar (UHV)). 

At normal conditions, equilibrium thermodynamics 
therefore predicts the observation of Ai4„_406n films^ 



We note that this result is not in contradiction with 
studies on the a-Al 2 O3(0001) or k-A1 2 C>3{001} surfaces, 
which predict Al termination! 19 i 32 In our study, the 
metallic substrate can take or give away charge via the 
interface, so that the polarity argument that rules out an 
O-terminated surface of a pure alumina slab cannot be 
applied to TiC/thin-film alumina. 

C. Trends in phase content, orientation, and 
preferred stacking 

It is clear that a detailed analysis of the trends in phase 
content, orientation, and preferred stacking of the alu- 
mina thin films must be based on the relaxed configura- 
tions. We find that a classification of thin- film candidate 
configurations that is based on the unrelaxed structures 
is dangerous. 

The stable and metastable alumina films are in gen- 
eral obtained from truncated TiC/K-Al2O3[001] interface 
configurations. This would be in agreement with the ex- 
perimental results that growth of K-AI2O3 is preferred 
over a-A^Os on clean TiC(lll) substrate^ However, 
truncated TiC/a-Al 2 03 and TiC/K-Al 2 O 3 [001] configu- 
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FIG. 4: Atomic structure of the stable A16O12 film. The color 
coding is: Dark gray = Ti, light gray = C, light = O, and black 
= Al. The top panel shows the projected side views along 
[100] and [010] including interlayer distances. The bottom 
panel shows the top views on the atomic layers [as defined in 
the top panel]. The Al coordination is OOO (O: octahedral, 
T: tetrahedral) . Note that the film is O terminated after 
relaxation. 



TABLE V: Stacking sequence and Al coordination of the un- 
relaxed alumina films with AUnOen stoichiometry and their 
relative total-energy differences E Te \ after relaxation. Nota- 
tion and grouping are the same as in Tab. Mil The ab ini- 
tio study and comparison permit a number of observations 
that are similar to those we made for the Al4n-406n films 
(Tab. IIII|I . However, the Al4„06n films are thermodynami- 
cally unstable, see Sec. IV Bl 



rations are (meta-)stable in the case of the thinner films 
with Ai4„_2 06n stoichiometry. 



Comparing configurations that are derived from struc- 
tures that differ only by a reflection about the yz-planc, 
we notice that, generally, the unrelaxed structure with 
AC stacking in the bottom two O layers yield more favor- 
able configurations than the structures with AB stack- 
ing. There are, however, exceptions; in particular the 
energetically most favorable Alio O is film and each one 
of the potentially metastable AlioOis and AI14O24 films 
posses an AB stacking in the bottom two O layers. We 
also note that, in general, the detailed stacking sequence 
of the energetically most favorable configurations varies 
strongly with the film thickness. 

In summary, although there are some general stabil- 
ity trends that can be inferred from the phase content, 
orientation, and stacking of the unrelaxed thin-film con- 
figurations, there are also several noticeable exceptions. 
In particular the Ai4„_2 06n films tend to break the rules. 
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FIG. 5: Atomic structure of the stable AI10O18 film (color 
coding as in Fig. [3}. The film is O terminated after relax- 
ation. The Al coordination is TjTj : T±00 (O: octahedral, 
T: tetrahedral, the arrows indicate the direction in which 
the tetrahedra point, different Al layers are separated by ':'). 
Note that tetrahedral coordination dominates and that both 
Al pairs in the bottom layer are tetrahedrally coordinated, 
with tetrahedra pointing in opposite directions (TfT|). 
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FIG. 6: Atomic structure of the two potentially metastable AlioOig films. Color coding and notation are as in Fig. [4] The 
Al coordinations are (a) OO : TjTjTj and (b) OO : T\TiO (O: octahedral, T: tetrahedral, the arrows indicate the direction in 
which the tetrahedra point, different Al layers are separated by ':'). 




VI. RESULTS II: ATOMIC STRUCTURE OF 
STABLE AND METASTABLE FILMS 



In this section we analyze in detail the atomic struc- 
ture of the relaxed alumina films with Ai4„_ 2 0g n and 
Ai4„_406r! stoichiometries that are found to be stable 
and potentially metastable in our implementation of the 
proposed method for structure search. While a num- 
ber of structural motifs that were included in the ini- 
tial candidate set are preserved, also novel motifs that 
strongly deviate from those in the bulk phase are iden- 
tified. Both preserved and novel motifs together, give 
insight into structural motifs of more accurate thin-film 
candidates. 



A. Atomic structure of the Ai4 n _2 06n films 

Two-O-layer thick films - Al§0\i. Figure [?] shows 
the atomic structure of the energetically most favorable 
AleOi2. It corresponds to a close-packed continuation of 
the TiC ABC substrate stacking, that is, the alumina 
stacking is Ab a bpb 1 C . All Al ions share the same atomic 
plane and are octahedrally (O) coordinated. The relaxed 
film is O terminated. Compared to TiC/O [O monolayer 
on TiC(lll)], the Ti-0 layer separation is drastically 
increased (+0.5 A). 

The potentially metastable AL3O12 structures possess 
almost the same structure as the energetically most fa- 
vorable one. They differ only by a slight displacement 



along the z direction of some of the Al ions. 

Three- O -layer thick films - AIiqOi%. Figure [5] shows 
the atomic structure of the energetically most favorable 
Alio O is film- Potentially metastable configurations are 
displayed in Fig. [51 In all cases, the Ti-0 layer separa- 
tion is shorter than in the energetically most favorable 
AlgOi2 film, but still considerably larger than in TiC/O 
(~ +0.3 A). We also note that in the stable AlioOis 
film, two of the six O ions in the bottom O layer are 
slightly lifted off from the TiC substrate. In the poten- 
tially metastable films, no O ion is lifted off. 

In all displayed films, the surface Al pairs have relaxed 
below the terminating O layer, so that the second Al layer 
consists of three Al pairs and the film is O terminated. 

The stacking of O layers is approximately described 
by ABA, AB(AC) bIidgQ1 and AC A for the energetically 
most favorable film and the two potentially metastable 
films respectively. The order is only approximate because 
a number of O ions are significantly distorted from ideal 
sites (as defined by the underlying TiC substrate). They 
are often located in bridge or cusp sites. This effect is 
most pronounced in the third O layer in the first poten- 
tially metastable film, which is entirely located in bridge 
sites. 

Similarly, the Al ions often deviate from ideal sites so 
that their description in terms of the bulk stacking labels 
becomes cumbersome. However, these distortions always 
occur pairwise, that is, Al pairs that are related by a 
bulk stacking label are dislocated symmetrically. The 
candidate structures generally preserve this symmetry of 
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FIG. 7: Atomic structure of the stable AI14O24 film (color 
coding and notation as in Fig. [4| . The film is O terminated 
after relaxation. The Al coordination is OT L : T T T| : OT4T4. 
(O: octahedral, T: tetrahedral, the arrows indicate the di- 
rection in which the tetrahedra point, different Al layers are 
separated by ':'). As in the most favorable AlioOis film, there 
is a layer in which both Al pairs are tetrahedrally coordinated, 
with tetrahedra pointing in opposite directions (T^Tf). Here, 
they are located in the second Al layer. 



the motifs of the bulk phases. 

The coordination of the Al ions is described as TfT± : 
T^OO, OO : T i T i T i and OO : OT t T h for the energeti- 
cally most favorable and the two metastable films respec- 
tively. Here and in the following O denotes octahedrally 
coordinated Al pairs, T tetrahedrally coordinated pairs. 
For tetrahedral coordination, T-\ means that the tetra- 
hedra point along the TiC[lll] direction, away from the 
interface, whereas Tj indicates that they point towards 
the interface. Different Al layers are separated by 

We notice that in all configurations, there is a large 
number of tetrahedrally coordinated Al ions (40-60%) 
and these can share the same atomic layer. In particular 
a larger number of tetrahedrally coordinated Al ions is 
favored. Furthermore, the energetically most favorable 
film contains tetrahedrally coordinated Al ions that share 
one layer and for which the tetrahedra point into opposite 
directions. 

Four-O-layer thick films - AI^Oia- Figures [7] and [5] 
show the atomic structure of the energetically most fa- 
vorable film and the two potentially metastable AI14O24 



films, respectively. The Ti-0 layer separations are com- 
parable to those in the energetically favorable Alio O is 
films. Also, in the stable AI14O24 film, two of the six 
O ions in the bottom O layer are slightly lifted off from 
the TiC substrate, whereas this is not observed in the 
potentially metastable films. In the most favorable film, 
the surface Al pairs have relaxed below the terminating 
O layer, so that the second Al layer consists of three Al 
pairs and the film is O terminated. However, both po- 
tentially metastable films are Al terminated even after 
relaxations. 

The stacking of O layers is more strongly distorted as 
in the case of AlioOis films, in particular in the most fa- 
vorable film. From the figures, we find the approximate 
O stacking sequences A(BC) bIidgB (A a AyCp)B [most fa- 
vorable film Fig. [7], ACBC [first potentially metastable, 
Fig. [51(a)], and ABAC [second potentially metastable 
film, Fig. [51(b)], where we use the the labeling (subscript) 
of the Al positions also for O ions and note that some of 
the ions are in fact dislocated from ideal sites. 

The Al ions are distorted correspondingly. This distor- 
tion is again pairwise and symmetrically for the energet- 
ically most favorbale and the first potentially metastable 
films. For the second potentially metastable film, this is 
not true. Both on the surface and in the first and sec- 
ond layer below the surface there are Al ions that have 
relaxed in a non-symmetric way. 

The coordination of the Al ions is given by OT^ : 
T X T T : OT i T i (most favorable), OT i : OO : OT i (first 
potentially metastable), and T\T^ : oTji^ : ooT^ (second 
potentially metastable). In the last sequence, the coor- 
dination of single ions that do not belong to a pair is 
denoted by small letters (t, o). Also, the coordination of 
the surface Al ions is not given for the two potentially 
metastable films. 

The result is similar to that for the AI14O24 films. In 
general, a large number of tetrahedrally coordinated Al 
ions is favored. In the most favorable film 70% of the Al 
ions are tetrahedrally coordinated. Furthermore, there 
is a layer with purely tetrahedrally Al ions and tetrahe- 
dra pointing into opposite directions (second Al layer). 
The first potentially metastable film possesses only 30% 
tetrahedrally coordinated Al ions. Inspection of the de- 
tailed stacking sequence, AbpayCciaapBbpa^Cba, identi- 
fies this structure as a partial K-AI2O3 configuration with 
an orientation TiC[Hl]/«[001]. 



B. Atomic structure of the Ai4„_406n films 

In general the potentially metastable Ai4 n _406ri films 
(if present) possess the same atomic structure as the the 
energetically most favorable film but rotated by 180° 
around the TiC[lll] direction [B <-> C,/3(2) <-► 7(3)]. 
They are therefore not discussed in the following. 

Two-O-layer thick films - Al^On- Figure [5] reports 
the calculated atomic structure of the energetically most 
favorable AI4O12 film. It is noticeable that the 0-0 sep- 
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FIG. 8: Atomic structure of the potentially metastable AI14O24 films. Color coding and notation are as in Fig. [4] The films are 
Al terminated even after relaxations. The Al coordinations are (without surface Al) (a) OTi : OO : TyO, (b) TyTy : Ty (tyo) : TyO 
(O: octahedral, T: tetrahedral, t single tetrahedral Al ion, o single octahedral Al ion, the arrows indicate the direction in which 
the tetrahedra point). 



aration in the alumina is relatively large, d _ ~ 2.6 A 
on average, and that the two Al pairs are not located 
between the O layers but are almost incorporated in the 
surface O layer, which leads to a large splitting of that 
layer. At the same time, the Ti-0 separation is compa- 
rably small and equals that in TiC/O. 

Hence, although predicted to be stable in a thermo- 
dynamical sense, structurally this TiC/Ai40i2 configu- 
ration separates into a TiC/O/ALjOg system, that is, a 
strongly bonded O monolayer on the TiC substrate with 
a thin alumina overlayer on top. 

The stacking of the O layers is AC, and the coordina- 
tion of the Al ions approx T^O. 

Three- O-layer thick films - Alg,Oi&. Figure [TOl shows 
the atomic structure of the energetically most favorable 
AlgOis fil m - The average 0-0 separation between the 
bottom two O layers is d _ ~ 2.5 A, which is slightly 
shorter than the one in the stable AI4O12 film. At the 
same time, the Ti-0 separation is increased to 1.15 A. 

The stacking of the O layers is approximately AC A. In 
the middle O layer, the two O ions that should be located 
in Cp are, however, dislocated to cusp sites. Furthermore 
the whole surface O layer is strongly distorted from ideal 
sites. The coordination sequence of the Al ions is T± : 



Only one of the original two Al pairs is left between the 
bottom two O layers after relaxation, the other pair mov- 
ing in between the top two O layers. Interestingly, one 
of the interfacial Ti atoms has left the Ti layer and re- 
laxed slightly in between the bottom two O layers. Again, 
structurally the TiC/AlsOis configuration appears as a 
partially decoupled TiC/0/AlgOi2 system. Here, how- 
ever, the Ti impurity above the bottom O layer may be 
a stabilizing factor. 

Four-O-layer thick films - AlyiO^- In Fig. 111! we 
show the atomic structure of the stable AI12O24 film. The 
O stacking is ACAB and hardly distorted. All Al ions 
have octahedral coordination. Thus, the present struc- 
ture mixes the O stacking of bulk K-AI2O3 with the Al 
coordination of bulk a-A^Os. 

The 0-0 separation d -o ~ 2.5 A is again very large 
and the Ti-0 separation is TiC/O like. Also, one of 
the original two Al pairs in the bottom Al layer has re- 
laxed upward through the middle Al layer and into the 
top layer. The other Al pair of the bottom layer is af- 
ter relaxation located only 0.1 A below the middle O 
layer. Consequently, also the TiC/Ali2024 configuration 
can again be considered as a decoupled, weakly binding 
TiC/0/Ali 2 0i 8 system. 
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FIG. 9: Atomic structure of the stable AI4O12 film (color cod- 
ing and notation as in Fig. [4| . Note the large 0-0 interlayer 
distance, which forms an almost empty region in between the 
two O layers. As a consequence, the TiC/alumina system sep- 
arates into TiC/O/alumina. The Al coordination is TjO (O: 
octahedral, T: tetrahedral, the arrow indicates the direction 
in which the tetrahedra point). 



VII. DISCUSSION 
A. Thin-film structure search method 

Although the relations between unrelaxed alumina 
structures and their energies after relaxations show some 
general trends, see Sec. IV C\ several important exceptions 
occur. In particular, these exceptions can result into the 
energetically most favorable structure. 

These exceptions illustrate a potential danger of apply- 
ing simple MC methods to the problem of finding the sta- 
ble thin-film oxide structures. An importance sampling 
of the thin-film configuration space based on a classifica- 
tion of the unrelaxed structures in terms of, for instance, 
alumina phase content, orientation, and/or O stacking 
may easily miss such exceptions. 

In our analysis of the relaxed atomic structure, we 
find that (apart from one possibly metastable AI12O24 
film) none of the stable and potentially metastable films 
shows a partial bulk alumina structure. Nevertheless in 
all structures the Al ions still obey parts of the bulk sym- 
metry. In particular, the positions of the Al pairs in the 
relaxed films are still pairwise related by the mapping 
given in Fig. [1] 

The set of thin-film candidates does have new struc- 
ture motifs different from those found in the bulk phases. 
The thin-film relaxation causes differences the coordina- 
tion of the Al pairs and/or in the stacking sequence of 
the O planes. For the Ai4„_206ri films we find that (i) 
layers with only tetrahedrally coordinated Al ions, T\T^ 
or Tj^Tj, are energetically favorable; (ii) layers with only 
octahedrally coordinated Al ions are present in both sta- 
ble and metastable films; (iii) Al layers with coordination 
TjTj are not present. For the A\± n -±OQ n films we find 
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FIG. 10: Atomic structure of the stable AlsOis film (color 
coding and notation as in Fig. [4| . Note that one of interfacial 
Ti atoms has relaxed to a position slightly above the bottom 
O layer. In the top view on OlAll in the lower panel, this Ti 
impurity is indicated by the large black ball. The Ti impurity 
may strengthen the TiC-alumina bond. The Al coordination 
is Tj, : OTyTi (O: octahedral, T: tetrahedral, the arrows 
indicate the direction in which the tetrahedra point). 



essentially only Al coordinations of the types OO and 
OTi- The TjTjO coordination of the Al layer directly be- 
low the surface O layer of the stable AlgOis film should 
be considered as a surface effect. The absence of OXf 
is consistent with the fact that the stable configurations 
all derive from TiC//c[001] sequences. However, the re- 
laxed configurations are still not conform with the bulk 
stacking. In particular, for the stable AI12O24 film, we 
observe only octahedrally coordinated Al ions in combi- 
nation with an ACAB stacking. 

The finding of new structure motifs, not explicitly in- 
cluded in the network of initial configurations, implies a 
significant strength. It shows that the proposed method 
is not restricted to a sorting of the original candidate 
structures in an energetic order, but that it is indeed ca- 
pable to predict energetically more favorable film geome- 
tries than what strictly constitutes symmetries in partial 
bulk structures. 

The identification of candidates for stable and 
metastable thin-film structures with some novel struc- 
tural motifs also suggests how the candidate space and 
the search could be broadened in a cost-efficient ap- 
proach. A broadening of the network of initial thin-film 
configurations can be made in the scope of structural el- 
ements. It would include motifs found in the bulk a- 
and K-AI2O3 but with a different weighting in choice of 
coordination for Al ions. In particular, the positions of 
the Al ions are still related pairwise, which restricts any 
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TABLE VI: Favorable structural motifs for relaxed films of different thicknesses. The label of the O layers indicate the 
approximate location of that layer with respect to the substrate surface (A — fee, B — hep, and C = top). The label of the 
Al pairs indicate their approximate coordination (O: octahedral, T: tetrahedral, the arrows indicate the direction in which the 
tetrahedra point). 
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FIG. 11: Atomic structure of the stable AI12O24 film (color 
coding and notation as in Fig. [3]). Note again the large in- 
terlayer distance between the bottom two O layers. The ab- 
sence of Al ions in between these O layers strongly indicates 
that the TiC/alumina system separates into weakly bound 
TiC/O/alumina. The Al coordination is O : OO : OOO, i.e., 
purely octahedral. 



broadening of the network of necessary initial thin-film 
configurations very significantly. We find indications for 
a significant increase in preference for tetrahedral coordi- 
nation of Al ions in the thin- film candidates. This motif 
is included in the initial thin-film network (which has 
many structures derived from K-AI2O3). However, a nat- 
ural further refinement of the present implementation of 
the proposed search method would be to include initial 



structures with a higher degree of tetrahedral coordina- 
tion of Al ions. It is possible to cast this broadening of 
the initial network into the framework of a genetic algo- 
rithm for identifying surface reconstructions ! 37 ' 38 ! 39 

We emphasize that a future, extended search for thin- 
film candidates is not expected to affect conclusions con- 
cerning thermodynamical stability of the various classes 
of ultra-thin alumina films. Since the slopes of T in 
Figs. [3] will remain unchanged, finding possibly energet- 
ically more favorable structures in the two relevant sto- 
ichiomctry classes will only resize the regions in which 
the different stoichiometrics are stabilized. To make 
the Ai4„_406n films generally unstable in comparison 
to Al4„_2 06n films, the truly stable Al4„_2 06n config- 
uration needs to gain at least ~ 10 eV compared to the 
stable Ai4„_406n found here. 



B. Note on the stability of CVD TiC/alumina 
wear-resistant coatings 

We emphasize that our results on the thermodynamical 
stability of thin-film alumina on TiC are critically based 
on the assumption of thermal equilibrium between the 
films and an O2 environment. The finding that the pre- 
ferred structure is described as a TiC/0/Al 4 („_ 1 )06(„-i) 
system (after relaxations) with a weak binding^ be- 



tween TiC/O and Al, 



4(n-l)06(n-l) 



is not in contradic- 



tion with the required exceptionally strong binding in a 
wear-resistant coating application. It rather shows that 
equilibrium between O2 and oxygen in the alumina dur- 
ing CVD growth of alumina is not reached. At environ- 
mental conditions relevant for CVD growth, as will be 
be discussed in a forthcoming paper ja! we find that it is 
instead the Al terminated films that arc stabilized. 



VIII. CONCLUSIONS 

We present a method to sample the configuration space 
of possible thin-film structures of complex oxides on a 
substrate. A well-defined network of initial configura- 
tions for promising thin-film candidates can be designed 
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from the oxide bulk structure. Ab initio calculations of 
relaxation deformations provide candidates for thin films 
as a function of stoichiometry and oxygen-layer thickness. 

The method has been illustrated for TiC/thin-film alu- 
mina, where experimental evidence^*^ can be used to 
reduce the network of initial thin films to contain struc- 
tural motifs defined by bulk a- and K-AI2O3. Based on 
this assumption, we have determined structural elements 
in and candidates for the energetically most favorable 
(stable or potentially metastable) TiC /thin-film alumina 
configurations for three thicknesses and three stoichiom- 
etry classes. 

Our method for ab initio search and study of thin-film 
structures has predictive power and provides detailed in- 
sight into the nature and atomic structure of thin-film 
alumina on TiC. The structures that are predicted by 
our method differ in their motifs heavily from motifs of 
the bulk structures, in particular in terms of the Al co- 
ordination. In principle, this warns that the present im- 
plementation of the search may not yet be complete and 
that we cannot make an authoritative prediction of the 
stable thin-film alumina structure; we can at present only 
identify key structural elements. More importantly, this 
finding of additional favorable motifs documents predic- 
tive power. It shows that the search method can identify 



candidates with a nature that is not explicitly included 
in the network of initial configurations. 

The different stoichiometry classes have been com- 
pared by means of Gibbs free energies. Assuming equi- 
librium with an O2 environment, we find that for the 
considered thicknesses of two, three, or four O layers 
(corresponding to n = 2, 3, or 4, respectively) the sta- 
ble films are either those with Ai4 n _406n stoichiometry 
(for medium to high O chemical potentials) or those with 
Al4„-206n stoichiometry (for very low O chemical poten- 
tials). The films with Al4„06n stoichiometry are never 
stabilized. 
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